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Running of the scalar spectral index in bouncing cosmologies 
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We calculate the running of the scalar index in the ekpyrotic and matter bounce cosmological 
scenarios, and find that it is typically negative for ekpyrotic models, while it is typically positive 
for realizations of the matter bounce where multiple fields are present. This can be compared to 
inflation, where the observationally preferred models typically predict a negative running. The 
magnitude of the running is expected to be between 10 -4 and up to 10~ 2 , leading in some cases to 
interesting expectations for near-future observations. 


I. INTRODUCTION 

Observations from Planck U! indicate that primor¬ 
dial perturbations were adiabatic, near-Gaussian, almost 
scale-invariant with a slight red tilt, and that tensor per¬ 
turbations were small. These observations have already 
greatly constrained the dynamics of the early universe. 

Our goal is to consider the running of the scalar spec¬ 
tral index as a further observational tool that can be used 
to differentiate between various cosmological scenarios, 
and in particular between inflation and two alternatives, 
the ekpyrotic universe and the matter bounce scenario. 
Even though the bounds obtained from observations al¬ 
ready do provide interesting constraints on these models, 
it will be helpful to have other observable quantities at 
our disposal to distinguish between these models, espe¬ 
cially since non-Gaussianities are already known to be 
small and that tensor modes have remained elusive. 

The primordial scalar perturbations are usefully char¬ 
acterised via their power spectrum 'Pn(k), which de¬ 
scribes the root mean square amplitude of the scalar cur¬ 
vature perturbations as a function of their wavenumber 
k. The spectral index of the perturbations is defined via 

din V-ji 


and it tells us how the amplitude of the perturbations 
is changing with scale. Here a denotes the scale factor 
of the universe, H = a/a the Hubble rate and we have 
defined A f as the number of e-folds left before the end of 
inflation, ekpyrosis or matter contraction, 

Af tot -Af = \n(a\H\). (2) 

For inflation, given that the Hubble rate H is nearly con¬ 
stant, one often defines the number of e-folds in terms of 
the variation of the scale factor alone. However, in bounc¬ 
ing cosmologies such an approximation is inappropriate, 
as the growth of \H\ is crucial. Note that in order for the 
flatness and horizon problems to be resolved, it is pre¬ 
cisely the quantity a\H\ that must grow by a factor e 60 
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or more before the hot big bang expansion phase. Ob¬ 
servations now indicate a clear deviation from a purely 
scale-invariant spectrum n s = 1, with Planck indicating 
a value [I], d] 

n s = 0.968 ± 0.006, (68% CL), (3) 

from the combined data from temperature fluctuations 
and lensing. Given that the spectrum has a slightly 
different amplitude at different scales, one may wonder 
whether the deviation from scale invariance also changes 
with scale. This is captured by the running a s of the 
spectral index, defined via 

dn s dn s 

= ITT = —T7’ ( 4 ) 

d n ^ at k=a\H\ 

Given that in all the models we are considering here 
the longest-wavelength modes get produced first, a s has 
the same sign as the time derivative of the spectral in¬ 
dex, dn s /dt. Such a running has not been observed yet, 
but Planck provides the following bound (again from the 
combined data from temperature fluctuations and lensing 

Hi) 

a s = -0.003 ± 0.007, (68% CL). (5) 

We are interested in how this measured bound compares 
with various theoretical models, and how future observa¬ 
tions may allow one to discriminate between models. We 
will begin by reviewing the predictions for inflationary 
models, and then contrast these with bouncing cosmolo¬ 
gies. Note that for the purposes of the present paper, 
we will consider the naive predictions of these cosmo¬ 
logical models and assume (perhaps optimistically given 
that these issues are not fully understood yet) that they 
are not significantly modified by eternal inflation or the 
properties of the bounce. 


II. REVIEW OF THE RUNNING OF THE 
SCALAR INDEX IN INFLATION 

The power spectrum of (single held) inhationary mod¬ 
els is given, up to a numerical factor, by 
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where H is the Hubble rate during inflation and e is the 
slow-roll parameter 
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here on, we will ignore corrections at sub-leading 
in e. Then, the scalar index n s is given by (see e.g. 


n s — 1 = — 
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while the running is 
d?i s 
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A very useful way of parameterising inflationary mod¬ 
els is obtained by combining the demand for a graceful 
exit (e = 1 when N = 0) with the assumption of a scaling 
symmetry Hi, 


1 

(N + 1)p ' 
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This parameterisation captures the qualitative behaviour 
of both chaotic (p = 1) and plateau ( p = 2) inflationary 
models, and it may be used to derive the predictions of in¬ 
flationary models with “minimal” tuning. With this pa¬ 
rameterisation, the spectral index and its running come 
out as 
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For chaotic models (pa 1), this leads to the typical val¬ 
ues n s ~ 0.95 and a s « — 8 x 10 -4 , assuming N = 60 for 
the modes of interest. For the plateau models favoured 
by current observations (p ~ 2), we find the typical val¬ 
ues n s ~ 0.97 and a s « —5 x 10~ 4 . In both cases (in fact 
in all cases captured by the above parameterisation) the 
running is negative in sign and of magnitude 0(1 /N 2 ) 
(note that in order to match the observed tilt in n s , p 
must lie between approximately 1 and 3). 

These results are simply presenting in a different fash¬ 
ion the well-known statement that, while inflation typ¬ 
ically predicts a negative running Q, a large negative 
running can be problematic for inflation [j]. That be¬ 
ing said, it is also known that if the inflaton potential 
contains modulations, this can increase the amplitude of 
the running (and even potentially give a positive run¬ 
ning) j8j. The drawback here is that this necessarily re¬ 
quires adding extra parameters to the theory, making it 
less predictive. Furthermore, the presence of sufficiently 
large modulations may make a graceful exit from infla¬ 
tion more difficult. 


III. THE EKPYROTIC UNIVERSE 


Instead of having an expansion phase with substan¬ 
tial negative pressure (inflation), the cosmological flat¬ 
ness problem may alternatively be solved by having a 
contracting phase with large positive pressure (ekpyro- 
sis). The high pressure suppresses both the homogenous 
curvature in the universe as well as small anisotropies, 
thus rendering the universe increasingly flat as it con¬ 
tracts @. A successful model then requires linking the 
ekpyrotic phase via a bounce to the current expanding 
phase of the universe — for implementations of smooth 
bounces see e.g. [Tol - ITtil and for a proof that the pertur¬ 
bations of interest pass through such bounces unaltered 
see [E [ 13 - 

In single-field ekpyrotic models, scalar perturbations 
are not amplified and thus the universe remains perfectly 
smooth, too smooth in fact! However, two-field models 
can produce scalar perturbations by amplifying isocur¬ 
vature/entropy modes first, and subsequently converting 
them into adiabatic curvature perturbations. Two types 
of such models have been considered to date, and we 
will analyse them in turn. For the first type, one in¬ 
troduces a potential which is unstable in the field space 
direction perpendicular to the ekpyrotic direction (which 
is the entropic direction) fl8l - l20l ]. This instability is re¬ 
sponsible for amplifying the entropic modes, but it also 
implies that the ekpyrotic phase now requires special ini¬ 
tial conditions (although in a cyclic context, this may be 
a feature rather than a problem Hi})- As detailed in j20| . 
the spectral index of the perturbations comes out as 


-| 2 q\ 

n s — 1 =-— . (13) 

e e 

Here e = — dlnff/dlna is defined in the same way as in 
inflation, with the difference that in the ekpyrotic case 
e is required to be larger than 3 (i.e., the potential is 
required to be steep). The calculation of the running 
then follows in a straightforward fashion, with the result 
that 
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For ekpyrosis, we can apply the same requirements of 
graceful exit (e = 3 when N = 0) and scaling symmetry 
to obtain the parameterisation Q 


e = 3(N + l) p . (15) 


This leads to the following expressions for the spectral 
index and the running: 
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Models fitting the observations well have p ~ 2, for which 
n s « 0.97. For these models, we may expect the running 
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to be given by a s ~ — 5 x 10 -4 . These predictions are 
in fact essentially identical to those of plateau models of 
inflation (also obtained for p ps 2). It is interesting to 
observe that in the present class of models all cases lead¬ 
ing to a red tilt (p ~ 1) also lead to a negative running 
of the spectrum. We should also remark that models 
of this ty pe can produce significant non-Gaussian correc¬ 
tions [221.1231]. which are in agreement with Planck bounds 
if the unstable potential is sy mmetric with respect to the 
ekpyrotic field direction [24J] . 

The second type of ekpyrotic model we will consider 
are models where the second scalar field does not acquire 
a potential (and remains stable), but where this scalar 
is non-minimally coupled to the ekpyrotic field <j>. The 
matter Lagrangians we have in mind are of the form ]25j ] 


£ = v=9 
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where the second scalar is denoted \ and it is kinetically 
coupled to the ekpyrotic field via a function fAs 
shown in (25| , for each ekpyrotic potential V (</) there ex¬ 
ists a coupling function f2(</) such that the perturbations 
in y obtain a specified spectrum, e.g. scale-invariant, but 
large departures from scale invariance are also possible. 
In fact, if we allow for arbitrary tuning, we may always 
find a function f 1 that yields a given spectral index and 
a given running. This is similar to the situation in in¬ 
flation, where, allowing arbitrary tuning of the potential, 
virtually any such combination may also be obtained. 
However, we are interested in models that are simple 
and that require as little fine-tuning as possible. For 
this reason, we will consider the case where the poten¬ 
tial can be approximated by the standard ekpyrotic form 
V{4>) = — Voe _c ^, but where we allow c to vary slowly 
over time so as to allow the ekpyrotic phase to have a 
graceful exit. Then for a coupling function H 2 = e~ b< ^ 
with b being a constant, the spectral index comes out as 


n s — 1 = —2 
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where we have assumed that e = c 2 /2 is large and we have 
dropped terms sub-leading in 1 /e. The situation we have 
in mind is where b and c are equal at the beginning of the 
ekpyrotic phase, and where c then gradually decreases 
according to e = c 2 /2 = 3(A7 + l) p . In this case, part 
of the observed red tilt can already be generated by the 
first term on the right hand side of m, implying that 
we would need to take p small, certainly p < 1. For this 
model, the running is given by 
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Again, the running is expected to be negative, and in 
fact it may be quite large in this class of models. For the 


extreme case where the entire red tilt is generated by the 
variation of e, we have p ~ 1 and thus a very large run¬ 
ning of a s ~ —10~ 2 , which is at a level that can be seen 
by near-future experiments. For the more general case 
where some of the red tilt is generated by the difference 
between b and c, we would expect p to be smaller and 
a s correspondingly also, but still at a significant level of 
a s ~ —10“ 3 for values of p within the range of validity 
of the approximations used here (since the spectral in¬ 
dex was derived using the approximation that e is large, 
we cannot take p smaller than approximately p = 1/10). 
Hence this type of model is special in that it predicts a 
much more significant running than the inflationary and 
ekpyrotic models discussed up to now, and this is some¬ 
thing to look out for in upcoming observations. 

Recently, new cosmological models have been proposed 
(anamorphosis [2?} and conflation [Hj) that combine in 
various ways elements of both inflation and ekpyrosis. 
Based on the arguments above we expect the running to 
also be negative for these models, while the amplitude of 
a s will depend on the specific realization. 


IV. THE MATTER BOUNCE SCENARIO 

As discovered by Wands [29}, a simple matter domi¬ 
nated contracting phase can also lead to scale-invariant 
scalar and tensor perturbations. This observation al¬ 
lows the construction of cosmological models in which 
the perturbations are generated in the presence of ordi¬ 
nary pressure-free matter [3d} . The horizon problem is 
automatically solved in matter bounce models since the 
big-bang singularity is resolved, and the flatness problem 
is alleviated: since the initial conditions are set far in the 
pre-bounce phase, it is possible to postulate a similar am¬ 
plitude of spatial curvature prior to the bounce point as 
is observed today (3lj. Note however that in a contract¬ 
ing FLRW space-time, a matter phase is not an attractor 
solution, and in particular the matter bounce is sensitive 
to the growth of anisotropies which in the absence of a 
fine-tuning of the initial conditions will come to dominate 
the high-curvature dynamics for generic solutions ]32,. [33} 
(although this last problem can be mitigated by adding 
an ekpyrotic phase after the matter-dominated era [H). 

In addition, and just as for ekpyrosis, the model must 
be completed by a description of a bounce. Realiza¬ 
tions of the matter bounce scenario typically fall into two 
categories depending on whether the bounce is caused 
by exotic matter fields that violate energy conditions 
13, [3a, |36[ or by modifications of the Einstein equations 
that are typically motivated by quantum gravity [37B45} . 
Note that the bounce can play an important role: for in¬ 
stance, bounces in loop quantum cosmology [39|,[45[ tend 
to suppress tensor fluctuations relative to scalar ones, 
which is a beneficial feature here as a contracting matter 
phase by itself tends to produce tensor fluctuations with 
an amplitude above present observational bounds j46} . 

Here, we will concentrate on the scalar fluctuations. 
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We will briefly review (and extend by allowing for a 
dynamical e) the calculation of their spectral index 
f39i Eli. [47, E|J in order to calculate the expected run¬ 
ning. The mode functions v = zlZ of the gauge-invariant 
curvature perturbations 1Z obey the Fourier space evolu¬ 
tion equation (for sound speed c s = 1) 

v" + (k 2 - v = 0 , (22) 


where a prime denotes a derivative with respect to the 
conformal time r, k is the comoving wavenumber and 
z = a(2e) 1 ^ 2 . It is easy to verify that, neglecting terms 
with two N derivatives on e (note that e"/2e contributes 
a term linear in e^), 
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where H = a'/a, dN = d ln(a) (or equivalently d M = (e— 
l)diV), and we have used the relation a"/a = V. 2 ( 2 — e) 
which follows from 0- 

An expression for H can be obtained from the defini¬ 
tion H' = H 2 ( 1 — e), which after successive integrations 
by parts (and neglecting terms of the order e^N and 
(e,Ar) 2 ) gives |2() 
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matter components with larger and larger equation of 
state successively come to dominate in a contracting uni¬ 
verse, we would expect the effective equation of state to 
grow during the matter phase, i.e. as M shrinks. Thus it 
is natural to expect wjj < 0 which would give a positive 
running, in contrast with inflation and ekpyrosis. 

A positive running is thus a typical signature of a mat¬ 
ter bounce model where multiple fields are present. How¬ 
ever, if there is only one field contributing to the dynam¬ 
ics, it is possible to obtain a negative running. One such 
example is given by models where it is a scalar field with 
an exponential potential that is acting as the pressure¬ 
less matter field (48|; the negative running arises in this 
case because, as mentioned above, the pressure-less solu¬ 
tion is not an attractor in a contracting space-time and 
in some cases the equation of state will shrink before 
growing again to reach the attractor which is kinetic- 
dominated [4E|. Nonetheless, in this paper we are inter¬ 
ested in “typical” realizations of the matter bounce and 
so long as there are a number of matter fields present in 
the space-time, one expects that the equation of state will 
generically increase in a contracting FLRW space-time. 

How large do we expect the running to be? The am¬ 
plitude of the running in the matter bounce varies signif¬ 
icantly from one model to another, but in any case the 
observational bound on a s gives the constraint 

|w,Af| < 10“ 3 , (68% CL). (30) 


Putting everything together (in the approximation that 
second order N derivatives are negligible), 
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Requiring the fluctuations to initially be in their vacuum 
state in the far past, the solution for the mode functions 
is given (up to an irrelevant phase) by 


v = {~ t) 7 ( 26 ) 

where H^ is a Hankel function with index 
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Rewriting this in terms of the effective equation of state 
of the matter fields to (related to e by w = |e — 1) and 
A f derivatives, we finally obtain 


n s — 1 = 12 oj — 9 wyvf, (28) 


and an immediate consequence is that the running can 
be approximated by 


a s = —12 w wf • 


(29) 


There is no simple way to estimate the quantitative be¬ 
haviour of the equation of state, unlike for the inflation¬ 
ary and ekpyrotic models above. However, given that 


One important consequence is that the main contribution 
to the red tilt of the scalar index, given in (128)) , must come 
from a negative w. Thus the spectral running provides 
both an interesting signature of the matter bounce, as 
well as a strong constraint for specific realizations of the 
matter bounce. 


V. DISCUSSION 

We have analyzed the predictions for the running of 
the scalar spectral index in ekpyrotic and matter bounce 
models, and contrasted these predictions with those ob¬ 
tained for simple inflationary models. The arguments 
given here are for “typical” realizations of the various 
cosmological scenarios. In that case we find that in all 
models the magnitude of spectral running is expected to 
be of the order of a few times 1/M 2 or about 1CU 3 for the 
modes of cosmological interest, and for non-minimally 
coupled ekpyrotic models and for the matter bounce the 
magnitude of the running may even be as large as 10 -2 . 
For inflationary and ekpyrotic models, the running is ex¬ 
pected to be negative, while for the matter bounce it is 
typically positive. Deviations from these predictions will 
require the models to be more complicated than the sim¬ 
ple parameterisations we have considered here (but might 
point to interesting missing ingredients), while an exper¬ 
imental confirmation of these predictions would provide 
a satisfying consistency test. 
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